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We present the first 3+1 dimensional simulations of non-Abelian plasma instabilities in gauge- 
covariant Boltzmann-Vlasov equations for the QCD gauge group SU(3) as well as for SU(4) and 
SU(5). The real-time evolution of instabilities for a plasma with stationary momentum-space 
anisotropy is studied using a hard-loop efFective theory that is discretized in the velocities of hard 
particles. We find that the numerically less expensive calculations using the group SU(2) essentially 
reproduce the nonperturbative dynamics of non-Abelian plasma instabilities with higher rank gauge 
groups provided the mass parameters of tlie corresponding hard-loop effective theories are the same. 
In particular we find very similar spectra for the turbulent cascade tliat forms in the strong-field 
regime, which is associated with an approximately linear growth of energy in collective fields. The 
magnitude of the linear growth however turns out to increase with the number of colors. 
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I. INTRODUCTION 

Non-Abclian plasma instabilities are, parametrically, 
the dominant coUective phenomenon in a weakly coupled 
quark-gluon plasma, and have been discussed as a possi- 
ble explanation for the extremely fast isotropization that 
is suggested by the success of hydrodynamical models 
of relativistic heavy-ion colHsions Such plasma in- 
stabilities are generalizations of thc so-called Weibcl or 
filamentary instabiUties in ordinary clcctromagnctic plas- 
mas . They are present already in coUisionless plasmas 
with any amount of momentum space anisotropy 0, 0] , 
and, indeed, they have been found to play an important 
role in the fast isotropization of electromagnetic plasmas 
Their non-AbeUan versions have been proposed to 
be of relevance for the quark-gluon plasma early on by 
Mrówczyñski and othcrs ¡6|-Í15¡ , and spccifically as expla- 
nation for thc fast apparent thermafization by Arnold ct 
al. [11, [13. 

NumericaUy, these instabifities have been studied in a 
discretized version of the hard-loop approximation [isl - 
[20| . which corrcsponds to gauge covariant Boltzmann- 
Vlasov equations describing thc dynamics of (soft) col- 
lcctivc ficlds in a weakly coupled plasma of hard parti- 
cles. The growth rate of plasma instabifities is parametri- 
caUy of the same order as plasma frequencies and screen- 
ing masses, and because of their exponential behavior 
the plasma instabifities dominate the coUective dynam- 
ics and incvitably lead to nonpcrturbatively large coUec- 
tivc ficlds. Eventuaüy thcy wiU givc risc to substantial 
backreactions on the momcntum distribution of thc hard 
particles, causing a breakdown of thc hard-loop approxi- 
mation ¡2l| - [23| coincident with the actual isotropization 
process. The hard-loop approximation aUows to study 
the early stage of this assumed scenario and thus its ba- 
sis. 

In thc carfiest stagcs of hcavy-ion coUisions it is in 



fact important to take into account thc cxpansion of the 
plasma, which modifies the growth from exponential in 
time to exponcntial in the square root of (proper) time, 
as indccd found in numerical simulations within the color 
glass condensate scheme [13, ÜHl as weU as in a general- 
ization of discretized hard loop simulations [26ll27|. In 
thosc simulations, an uncomfortable delay of the onset 
of growth has been observed, which however has recently 
been shown to largely disappear upon consideration of 
more general initial conditions [l^. While the density 
and fifctimc of thc plasma cstimated for heavy-ion col- 
fisions at RHIC may bc too low to give an important 
role to nonabelian plasma instabilities there, the higher 
values cxpected for LHC hcavy-ion coUisions may be suf- 
ficient for nonabeUan plasma instabifities to become the 
dominant phcnomcnon in a lcss strongly coupled envi- 
ronment. 

By means of real-time lattice simulations for stationary 
anisotropic plasmas, it has however been found that in 
contrast to effectively 1-1-1 dimensional situations where 
only the most unstable modcs arc considercd, thc cxpo- 
nential growth of non-Abcfian plasma instabifities is fim- 
itcd in 34-1 dimcnsions by non-Abcfian seU-interactions. 
At a certain magnitude of the nonabefian fields, which 
depends on the degree of anisotropy [Í^, , the expo- 
nential growth ceases and turns into a finear growth of 
tfie energy densities of the soft ficlds. In that regime, 
a turbulent cascade of energy is observed, with a quasi- 
steady-statc power-law distribution oc k^'^ and spec- 
tral index j/ 2, which transports thc energy fed into 
low-lying modes by Wcibcl instabifitics to stable highcr- 
momentum (plasmon) modes through non-Abefian seU- 
interactions of gluon fields [3l|. This cascade forms at 
momentum scales that are parametricaUy separated from 
those of the hard particles, making it possible to study 
this phenomenon seU-consistently within the hard-loop 
approximation. On the other hand, classical-statistical 
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simulations in SU(2) gauge theory [32, Is^, whcre there 
is no such separation of scales, have reported a late-time 
behavior indicative of a lower index of = 4/3 (as in 
turbulence with constant transport of particle number) 
while showing an early behavior quahtatively similar to 
those of Chromo-Weibel plasma instabihtics. The index 
V = 2 found in the hard-loop simulations was argued ¡33 | 
to fit to an encrgy cascadc carried by particles in the cas- 
cade scattering off of nonperturbatively large background 
fields. However, this is not consistent with a steady-state 
particlc cascade [s^ (particle number of cascade particles 
is unchanged in the assumed process) so that either other 
processes are equahy important or the resuhing cascade 
is not steady-state. (See also Refs. [Ü, [H, Isa - li^ for 
discussions of instabihties and turbulence in QCD.) 

In this paper we confirm the resuhs [l^ [13, l29l - [3ll | ob- 
tained for stationary anisotropic plasmas in the hard-loop 
framcwork for gaugc group SU(2), and consider, for the 
first time, 3-(-l-dimensional hard-loop simulations for thc 
QCD gauge group SU(3), as weh as for SU(4) and SU(5), 
in order to quantify the dependence on the number of col- 
ors. So far, in the hard-loop effective theory, SU(3) cal- 
culations have only been performed for effectively 1-1-1 
dimensional situations [l^. (In classical-statistical 3-1-1- 
dimensional lattice gauge theory the gauge group SU(3) 
has becn studied in Ref. |43|.) 



II. SETUP 

In order to study nonabehan plasma instabihties in a 
wcakly coupled quark-gluon plasma, which may for the 
first time become physicaUy relevant at the higher tem- 
peratures and densities reached in heavy-ion cohisions at 
the LHC, we consider thc extreme hmit of an uhrarel- 
ativistic coüisionless plasma, where the main dynamics 
takes place at scales parametricaUy soft compared to the 
hard scale |p| = of the plasma constituents. In an 
isotropic plasma, the scale gjpl, where g is the gauge cou- 
phng, determines the scale of the Debye screening mass, 
of the plasma frequency, and of Landau damping, but 
in an anisotropic plasma, the dominant cohective phe- 
nomenon at this scale, which is larger than the scale of 
coUisions, turns out to be plasma instabihties. 

The effective field theory relevant for the coUective phe- 
nomena at this largest of the soft scales is given by gauge- 
covariant cohisionless BoUzmann-Vlasov equations ¡44| 
or "hard-loop" effective theory as long as the soft gauge 
ficlds obey ^ |p|/5 so that a backrcaction on the hard 
particles can stiU be ignored. This aUows one to study 
in detail the first stage of plasma instabihties, which, as 
briefly reviewed in the introduction, is highly nontrivial 
in nonabehan gauge theories. 

The corresponding effective action, which is nonlocal 
and nonhnear ¡45l . |46¡ , can be made local at the expense 
of introducing auxihary fields [47| in the adjoint represen- 
tation, W^(a;;v), for each spatial unit vector appcaring 
in thc velocity ti^ = p^/\\i\ = (l,v) of a hard (uUra- 



relativistic) particle with momentum p^ . The W ficlds 
encode the fluctuations of the distribution function of 
colored hard particles. In terms of these, the induced 
current j in the nonabehan MaxweU equations 



¿?^(yi)F^'' = r, 



(1) 



is given by 



^ J (27r)3 2|p|^ dpf^ ^ ' ^ ' 



and the nonabchan Boftzmann-Vlasov equation, in which 
the scale of the uUrarelativistic hard particles drops out, 
reduces to 



[v ■ D{A)]Wp{x;^) = Fp^{A)v\ 



(3) 



with Dfj^ = df¿ — ig[Af¿, ■]. (Our metric convention is (H — 
— )•) 

We shaU only consider a background distribution 
ñmction /(p) of hard particles with one direction of 
anisotropy, obtaincd by deforming an isotropic distribu- 
tion according to 



/(p)«/íso(p'+?pD 



(4) 



which leads to 
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- [a{yr)W"{x; v) + 6(v)VK^(x; v)] v^" 
W{x;x)v>', (5) 



and 



a(v) 



(1+^1^2)2' 



b{v) = ^Vza{-v), 



(6) 



with the mass paramcter m proportional to g and the 
scale of the momenta of hard particles. 

In the numerical treatment we discretize the 3- 
dimensional configuration space by cubic lattices, and the 
unit sphere by a set of unit vectors v pointing to 

= -1 + (2¿ - l)/iV„ z = l...iV„ 
(^, = ^TTj/iV^, j = l...N^ (7) 

in cyUndrical coordinates. At each lattice site we thus in- 
troduce A/W = J^z xAL variables Wv{x) whose dynamics 
is approximated by [13, HO] 

[v ■ D{A)]Wv = (avJ^°'^ + bvF'-^'^v^ (8) 

D,{A)F'^'' = f = ^Y.^'^Wv (9) 

V 

A different possibihty of discretization was used in 
Ref. [l^, where the auxihary fields W^{x;v) are ex- 
panded in spherical harmonics Yi„i{-v) and truncated at 

SOme Itn.n^r- 
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In thc foUowing we shall considcr a momentum space 
distribution of hard modes which is oblate, by choosing 
^ = 10 as beforc in Refs. [3 H^l- In this situation, a 
prolate region in momentum space at soft scales, pinched 
at the origin and involving longitudinal momentum < 
\kz\ < M = 1.072m, contains unstable modes. The mode 
with largest growth rate 7* = 0.12m is found at = 
and fcz = fc* = 0.398m. The mass paramcter m dcfined 
above is related to the asymptotic mass of propagating 
transverse gluons by moo = 0.447to. 

The details of the lattice discretization of the above 
nonabehan gauge-covariant Bohzmann-Vlasov equations 
can be found in the appendix of Ref. [l^. The only 
place where those formulae need to be modified to ahow 
for general is in the expression for the energy density. 
For arbitrary N^, the chromo-magnetic part of the energy 
density in Eq. (B18) of Rcf. rcads 

^^ = ^E(i-¿t'^t^ü(*'-))- (10) 

In the numerical simulations of the above equations 
we have used paraüehzation on computer clusters, spht- 
ting up thc simulation domain so as to fit into the main 
memories of the individual nodes, and considered total 
spatial lattices ranging from 32^^ to 128^ and numbers of 
W fields TVw = A/; X A/'y = 20 X 16 = 320 and higher 
at each site, in addition to the Ng = N^ — 1 gluon fields 
(fink variables) and their conjugate momenta. 

Thc physical size of the lattice is detcrmincd by the 
parametcr rn^ appearing as prcfactor in the induced 
currcnt, Eq. ([S]). In terms of the asymptotic thermal 
mass of gluons. thc latticc spacings that we have con- 
sidered vary from amoo = 0.5656 to 0.1414. For the 
specific anisotropy that wc are considering, this corre- 
sponds to a/(A*/2) = 0.16 . . .0.04 in terms of the wave 
length of maximal growth A,, or, with respect to the 
minimal wave length A^ for which unstable modes exist: 
a/(A^/2) = 0.432... 0.108. 

In order to providc sced fields, we have initiahzed with 
vanishing gauge fields and 

(>V:(0,x)W^(0,j/)) = S'^'Siya^ (11) 

where 6^ y dcnotes 3 Kroneckcr dcltas, subtracting 
'^^W"' /Aíw'^ from each W"' in order to obey Gauss's 
law. 

III. NUMERICAL RESULTS 

A. Growth of energy densities in 
chromo-electromagnetic flelds 

In the foUowing we shah concentrate on the case of 
small initial ficlds (small a) so that there is first a phase 
of esscntiaUy Abchan evohition with exponential growth 
due to the Wcibcl instabihty, which ends whcn thc gaugc 




FIG. 1: (Color online.) Comparison of average energy densi- 
ties £ for SU(2) (light) and SU(3) (dark) on logarithmic (up- 
per panel) and linear (lower panel) scale in 3-1-1 dimensional 
simulations for anisotropy parameter ^ = 10 on a 64'^ lattice 
with ajTioo = 0.5656 (i.e., a = 0.16(A*/2), L = 64a = 13.8A^) 
and TVw = 320 (= Nz x N^ = 20 x 16). The initialization 
of the W fields is chosen with parameter a = 0.1 resulting in 
equal initial field amplitudes per gluon degree of freedom for 
both gauge groups. Shown are tlie total field energy densities 
as well as the contributions from transverse and longitudi- 
nal chromo-electric and chromo-magnetic fields. In the lower 
panel, the individual runs are shown by light curves. 

fields have grown to nonperturbatively large ampfitudes 
so that non-Abefian self interactions become important. 
We keep fixed the anisotropy parameter ^ = 10 in aU of 
the foUowing comparisons. 

For initiaUy smaU fields, Figurc [T] shows a comparison 
between SU(2) and SU(3) simulations. Thc simulations 
have been initiafized with random fiuctuations in thc W 
field with strength a = 0.1 for both gauge groups. De- 
picted is the total energy density and its contributions 
from chromo-electric and chromo-magnetic fields, each 
decomposed into parts transverse and longitudinal with 
rcspcct to the direction of anisotropy. The SU(2) result 
has bccn obtained from averaging over 6 runs with diffcr- 
ent initial conditions, while for SU(3) 3 runs havc bcen 
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FIG. 2: (Color online.) Same as Fig. [T]for the energy density 
divided by the number of gluons Ng. 



used. Because there are 8 gluons in SU(3) comparcd to 3 
gluons in SU(2) that are initiahzed with the same avcragc 
gauge field amphtude, the initial total energy density of 
SU(3) is larger than the SU(2) energy density by a factor 
8/3. Remarkably, at the end of the exponentiaUy growing 
phase around í ~ 60 when the approximately hnear 
growth regime sets in, the total cncrgy densitics approx- 
imatcly agree (after averaging ovcr individual runs — the 
saturation energy of individual runs varies substantiahy, 
in particular for smaUer lattice size). We have checkcd in 
comparison runs restricted to an Abehan U(l) subgroup 
that this agreement is not caused by the lattice spac- 
ing chosen and that we are safely below the compactness 
bound. The plot in hnear scale shows that the SU(3) re- 
suh subscqucntly grows somewhat faster than thc SU(2) 
result. 

In Fig. [2] the same results are compared by dividing 
by the numbcr of ghions Ng. Now thc SU(2) and SU(3) 
curves agree in the exponentiaUy growing phase, with the 
SU(3) result starting to lag behind the SU(2) resuh from 
moo í ~ 35 onwards. In the hnear growth regime, it is 
now thc SU(2) curve that grows faster. Thercforc thc 
dcpendcncc of the hnear growth rate in the non-Abchan 



FIG. 3: (Color online.) Comparison of average total field 
energy densities £ for SU(2) through SU(5) on logarithmic 
(upper panel) and linear (lower panel) scale in 3+1 dimen- 
sional simulations for anisotropy parameter ^ = 10 on a 
64=* lattice with A^w = 320 (= A^z x A^^ = 20 x 16). 
The initialization of the W fields is chosen with parameter 

0"SU{2) = fSUCS) = 0-SU{4) = f"SU(5) = 0.1. 

regime on N^ sccms to bc a factor betwecn 1 and Ng. 

In order to study the systematics of thc scahng with 
Nc in the non-Abehan regime, we extend the gauge group 
SU(A''c) to larger values of Nc. Figure[3]shows a compari- 
son for gauge groups SU(2) through SU(5). This confirms 
thc obscrvation of Fig.[T]that thc cncrgy dcnsitics stop to 
grow at approximately thc same valuc for different gauge 
groups, and also that the energy densities of the larger 
gaugc groups grow faster in the regime of hnear growth. 

In Fig.|3]we plot thc ratio of the time dcrivatives of the 
energy densities of Fig. [3] for thc various gaugc groups, 
ah with same lattice parameters and initial gaugc ficld 
strcngths. Sincc thc discrete derivatives wigglc strongly, 
the energy densities have been smoothed by averaging 
over a rectangular window function of size mooAí = 2. 
These resuhs indicatc that thc slope of the energy density 
in thc hncar growth rcgimc scalcs approximatcly propor- 
tional to thc numbcr of colors, N^. 



5 




SU(3)/SU(2) 
SU(4)/SU(2) 
SU(5)/SU(2) 




SU(3)/SU(2) 
SU(5)/SU(2) 



FIG. 4: (Color online.) The ratio of the time derivative of 
the energy densities of Fig. [3]for gauge groups SU(3), SU(4), 
and SU(5), divided by the one of SU(2). For large times, this 
is seen to scale approximately with Nc- Horizontal Unes at 
3/2, 4/2, and 5/2 mark the ratios for scaUng by the number 
of colors. 



FIG. 5: (Color onUne.) The ratio of the time derivative of 
the energy densities for various lattice parameters. (a) 32"^, 
am^ = 0.2828, a = 0.0707; (b) 32^ amoo = 0.2828, a = 3.46; 
(c) 64^ amoo = 0.1414, a = 0.1. For the SU(5)/SU(2) slope 
ratio, the same parameters as for (b) have been used. 



In Fig.[5]we show the slope ratio of gauge groups SU(3) 
and SU(2) for finer lattices (higher ultraviolet cutoff) at 
equal (but now smaller) physical volume size. For these 
curves, a rectangular window function of size TOoo Aí — 10 
has been used. The data summarize averages over (a) 20, 
(b) 8, (c) 3 runs with diffcrcnt random initial conditions 
for SU(2), (a) 3, (b) 4, (c) 1 runs for SU(3), and one run 
for SU(5). In these runs, which have less resolution in thc 
infrared, we stih find ratios that are consistently larger 
than 1, but more scattered about the ratio of the number 
of colors. A scahng by number of colors is also supported 
by the slope ratio of SU(5) and SU(2) in Fig. É We 
are therefore led to the conjecturc that thc finear growth 
rate in the non-Abefian regimc is proportional to Nc, but 
have to concede that our numcrical vcrification involves 
sizeable uncertainties. 



B. Spectra 

In order to get a gfimpse of tfic underlying dynamics, 
fofiowing Ref. [3l| we considcr distribution functions of 
modes 



fAÍk) = 
ÍEÍk) = 



1 



NdoikV 



(A2(fc)), 
[-E^k)}, 



(12) 



where V is the total spatial volume and A'^dof = 2A^g ac- 
counts for two transverse polarization states and Ng = 
N^ — 1 adjoint color states of the SU(iVc) gauge theory. 
The spectra are obtained from a Fourier transformation 
of the A and E ficlds in tfie lattice Coulomb gauge ¡48| . 
which fixes the residual gauge freedom of the temporal 



axial gauge by minimizing unphysical high-momcntum 
noise within 3-dimcnsional time sfices. Tfiis turns out to 
be essential for our resufis for the spectra which would 
otherwise show much more power in the ufiraviolet. Sta- 
ble plasma modes are expected to contribute equaUy to 
/a and /e, whereas unstable modes are predominantly 
magnetic, leading to /a > Íe in thc corresponding mo- 
mcntum range. 

Figure [6] shows the corresponding spectra, normafized 
per gluonic dcgree of freedom. The results are averaged 
over lattice vectors k in 12 equally sized bins on a loga- 
rithmic scale along |k^|. In addition, SU(2) curvcs rep- 
resent averages ovcr 6 runs, while SU(3) curves represent 
averages over 4 runs. 

Thc upper pancl of Fig. [Ü] compares the magnetic 
(/^(fc)) and electric ifsik)) distribution functions. The 
lowest curves correspond to the initial timc. Around 
k/moo ~ 0.9 there is a conspicuous cxponential growth 
of modes associated with Weibel instabilities. In order 
to ease the comparison of /a and /e results, the time 
mooAí « 14 is singled out in this plot by thick fines. As 
can be seen also in Fig. [U the magnetic contributions 
dominate over electric contributions during the cxponen- 
tial growth phase. Whcn non-Abelian sefi-interaction of 
the gluonic fields sets in, tfie exponential growth ceases, 
and only the higher modes fc/moo ^ 2 continue to grow, 
albeit more slowly. 

A comparison between SU(2) and SU(3) is shown in 
the lower panel of Fig. [üj As expected, the behavior of 
the exponentially growing modes agrees well between the 
SU(2) and SU(3) calculations since in that regime non- 
Abefian seff interactions are for the most part smafi. Note 
that the definitions for fAÍk) and fsik) include a 
scafing by the number of degrees of freedom, which is 
proportional to Ng. This plot therefore corresponds to 
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FIG. 6: (Color online.) Spectra corresponding to Fig. [T]for 
times < rrioat < 150. The distance between the lines is 
mooAí 6. The upper panel shows the SU(2) spectra for 
the magnetic (/^(fc)) and electric {fE{k)) distribution func- 
tions. The curves at time mooí ~ 23 are plotted with thick 
lines. The lower panel shows the electric (/^(fc)) distribution 
functions for SU(2) vs. SU(3). 



the energy density divided by Ng as depicted in Fig. [2] 
The evolution of thc SU(2) and SU(3) curves only deviate 
from each other when non-hncaritics sct in. 

In order to study the non-Abehan regime with approxi- 
mately hnear growth of energies more carefuUy, Fig.[7]de- 
picts the spectra multiphed by fc^ and at late-time times 
80 < moot < 150. The distribution functions have been 
averaged over 32 bins along |k^| in these plots so that 
the large k rcgion is better resolvcd. Thc slow growth 
at large momcnta k/rrioo ^ 2 corresponds to thc linear 
growth rcgimc of Fig. [1] On larger scales k/moo ^ 2 
we have ¡A^k) — ¡E^k) as expected from stable modes 
in the perturbative regime. In an intermediate range 
0.5 < k/moa ^ 2, magnetic fields dominate over electric 
ones fA^k) > fsik), which is caused by the predomi- 
nantly magnetic plasma instabilities. 

The straight hnes in the upper panel of Fig. [7] cor- 
respond to simplc power laws f{k) ~ k~^ with v = 2 
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FIG. 7: (Color onhne.) The power spectrum for SU(2) of the 
Coulomb gauge distribution f{k) at late times 80 < rrioot < 
150. The distance between the lines is mooAí « 11. The up- 
per panel shows a comparison for SU(2) for fs^k) and fA{k). 
The straight central horizontal red line indicates a power law 
spectrum / ~ fc^" with p = 2, while the dotted lines cor- 
respond to v = 1.8 and 2.2. The lower panel compares the 
spectra fe^k) for SU(2), SU(3), SU(4), and SU(5). 



(straight hnc) and i/ = 1.8, and 2.2 (dottcd hncs). A dis- 
tinct power law bchavior k^^ with v k, 2 over a large 
range of momenta can be most clearly extracted from 
the electric distribution function fs^k), which does not 
show a hump caused by thc unstablc magnctic modcs, 
and if onc cxcludes the higher momcntum modcs closc to 
thc latticc cutoff and also the largcst timcs, whcrc cffects 
from thc latticc cutoff may be felt. As error bar for i/ w 2 
we infcr 1.8 < < 2.4. 

The lower panel of Fig. [7] compares the electric ficld 
spectra fs^k) at various gaugc groups SU(2) through 
SU(5). Thc spcctra of the various gaugc groups appear 
to bc similar, apart from an ovcrah scahng factor. 
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IV. CONCLUSIONS 

We have studied for the first time 3+1 dimensional sim- 
ulations of non-Abehan plasma instabihties in the hard- 
loop framework for the gauge group SU(3) that is phys- 
icany relevant for the physics of heavy-ion coUisions. In 
order to study the dependence on the number of colors, 
we have also considcred gaugc groups SU(4) and SU(5). 
We have found that smah seed fields which correspond 
to the same amount of initial energy density per gluon 
degree of freedom lead to comparable total energy den- 
sity when exponential growth of plasma instabihties is 
stoppcd by nonabehan self-intcractions and a phase of 
approximatcly lincar growth bcgins. Wc havc confirmcd 
that in aü nonabehan gauge groups considered a power- 
law spectrum /(fc) ^ with f 2 develops for the 
higher-momentum modes, corresponding to a cascade of 
energy towards the ultraviolet. The growth rate of en- 



ergy dcnsities in this phase was found to scalc roughly 
proportional to the number of colors. 
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